Coupled spin and valley physics in monolayers of M0S2 and other group-VI dichalcogenides 
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We show that inversion symmetry breaking together with spin-orbit couphng leads to coupled spin and valley 
physics in monolayers of M0S2 and other group-VI dichalcogenides, making possible controls of spin and 
valley in these 2D materials. The spin-valley coupling at the valence band edges suppresses spin and valley 
relaxation, as flip of each index alone is forbidden by the valley contrasting spin splitting. Valley Hall and 
spin Hall effects coexist in both electron-doped and hole-doped systems. Optical interband transitions have 
frequency-dependent polarization selection rules which allow selective photoexcitation of carriers with various 
combination of valley and spin indices. Photo-induced spin Hall and valley Hall effects can generate long lived 
spin and valley accumulations on sample boundaries. The physics discussed here provides a route towards 
the integration of valleytronics and spintronics in multi-valley materials with strong spin-orbit coupling and 
inversion symmetry breaking. 
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Since the celebrated discovery of graphene [1-3], there has 
been a growing interest in atomically thin two-dimensional 
(2D) crystals for potential applications in next-generation 
nano-electronic devices [4, 5]. Layered transition-metal 
dichalcogenides represent another class of materials that can 
be shaped into monolayers [4], which display distinct physical 
properties from their bulk counterpart [6-9]. Recent exper- 
iments have demonstrated that M0S2, a prototypical group- 
VI dichalcogenide, crossovers from an indirect-gap semi- 
conductor at multilayers to a direct band-gap one at mono- 
layer [6, 7]. The direct band-gap is in the visible frequency 
range, most favorable for optoelectronic applications. Mono- 
layer M0S2 transistor was also realized, demonstrating a 
room-temperature mobility over 200 cm^/(V-s) [8]. 

In monolayer M0S2, the conduction and valence band 
edges are located at the corners {K points) of the 2D hexag- 
onal Brillouin zone [10-12]. Similar to graphene, the two 
inequivalent valleys constitute a binary index for low energy 
carriers. Because of the large valley separation in momentum 
space, the valley index is expected to be robust against scat- 
tering by smooth deformations and long wavelength phonons. 
The use of valley index as a potential information carrier was 
first suggested in the studies of conventional semiconductors 
such as AlAs and Si [13]. With the emergence of graphene, 
the concept of valleytronics based on manipulating the valley 
index has attracted great interests [14-18]. 

M0S2 monolayers have two important distinctions from 
graphene. First, inversion symmetry is explicitly broken in 
monolayer M0S2, which can give rise to the valley Hall effect 
where carriers in different valleys flow to opposite transverse 
edges when an in-plane electric field is applied [15]. Inversion 
symmetry breaking can also lead to valley-dependent optical 
selection rules for inter-band transitions at K points [16]. Sec- 
ond, M0S2 has a strong spin-orbit coupling (SOC) originated 
from the d-orbitals of the heavy metal atoms [12], and can 



be an interesting platform to explore spin physics and spin- 
tronics applications absent in graphene due to its vanishing 
SOC [19, 20]. 

In this Letter, we show that inversion symmetry break- 
ing together with strong SOC lead to coupled spin and val- 
ley physics in monolayer M0S2 and other group-VI dichalco- 
genides, making possible spin and valley control in these 2D 
materials. We find the conduction and valence band edges 
near K points are well described by massive Dirac fermions 
with strong valley-spin coupling in the valence band, which 
has several important consequences. First, the valley Hall ef- 
fect is accompanied by a spin Hall effect in both electron- 
doped and hole-doped systems [21-24]. Second, spin and val- 
ley relaxation are suppressed at the valence band edges as flip 
of each index alone is forbidden by the valley-contrasting spin 
splitting (~ 0. 1 - 0.5 eV) caused by inversion symmetry break- 
ing. Third, the valley-dependent optical selection rule also be- 
comes spin-dependent, and carriers with various combination 
of valley and spin indices can be selectively excited by optical 
fields of different circular polarizations and frequencies. We 
predict photo-induced charge Hall, spin Hall and valley Hall 
effects. The latter two phenomena can be used to generate 
long lived spin and valley accumulations on sample bound- 
aries. The physics discussed here provides a route towards 
the integration of valleytronics and spintronics in multi- valley 
materials with strong spin-orbit coupling and inversion sym- 
metry breaking. 

The physics in monolayers is essentially the same for 
group-VI dichalcogenides MX2 (M = Mo, W, X = S, Se), de- 
scribed below using M0S2 as an example. Structurally, M0S2 
can be regarded as strongly bonded 2D S-Mo-S layers that are 
loosely coupled to one another by Van der Waals interactions. 
Within each layer, the Mo and S atoms form 2D hexagonal 
lattices, with the Mo atom being coordinated by the six neigh- 
boring S atoms in a trigonal prismatic geometry (Fig. la-b). In 
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FIG. 1. (color online), (a) The unit cell of bulk 2i7-MoS2, which has 
the inversion center located in the middle plane. It contains two unit 
cells of M0S2 monolayers, which lacks an inversion center, (b) Top 
view of the M0S2 monolayer. Ri are the vectors connecting nearest 
Mo atoms, (c) Schematic drawing of the band structure at the band 
edges located at the K points. 



its bulk form, M0S2 has the 2H stacking order with the space 
group Dq^, which is inversion symmetric. When it is thinned 
down to a monolayer, the crystal symmetry reduces to Dl^, 
and inversion symmetry is explicitly broken: taking the Mo 
atom as the inversion center, an S atom will be mapped onto 
an empty location. As a consequence, the effects we predict 
here are expected only in thin films with odd number of lay- 
ers, since inversion symmetry is preserved in films with even 
number of layers. 

We start by constructing a minimal band model on the 
basis of general symmetry consideration. The band struc- 
ture of M0S2, to a first approximation, consists of partially 
filled Mo (i-bands lying between Mo-S s-p bonding and anti- 
bonding bands [25]. The trigonal prismatic coordination of 
the Mo atom splits its (i-orbitals into three groups: Ai{dz2), 



E{dxy^ dx2-y2) and E'{dxz^ dyz). In the monolayer limit, the 
reflection symmetry in the z direction permits hybridization 
only between Ai and E orbitals, which opens a band gap at 
the K and —K points [25], schematically shown in Fig. Ic. 
The group of the wave vector at the band edges {K) is C^h 
and the symmetry adapted basis functions are 
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where the subscript c{v) indicates conduction (valence) band, 
and r = ±1 is the valley index. The valence-band wave func- 
tions at the two valleys, |^+) and |^~), are related by time- 
reversal operation. To first order in k, the C2,h symmetry dic- 
tates that the two-band k • p Hamiltonian has the form 



^0 = atirkxCFx + kyCFy) + — 



(2) 



TABLE I. Fitting result from first-principles band structure calcu- 
lations. The monolayer is relaxed. The sizes of spin splitting 2 A 
at valence band edge were previously reported in the first principle 
studies [12]. The unit is A for a, and eV for t, A and A. Qi (Q2) is 
the Berry curvature in unit of A^, evaluated at — point for the spin 
up (down) conduction band. 
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M0S2 


3.193 


1.66 


1.10 


0.15 


9.88 


8.26 


WS2 


3.197 


1.79 


1.37 


0.43 


15.51 


9.57 


MoSe2 


3.313 


1.47 


0.94 


0.18 


10.23 


7.96 


WSe2 


3.310 


1.60 


1.19 


0.46 


16.81 


9.39 



where cr denotes the Pauli matrices for the two basis functions, 
a is the lattice constant, t the effective hopping integral, and 
A the energy gap. These parameters are obtained by fitting 
to first-principles band structure calculations and are listed in 
Table. I for the four group- VI dichalcogenides [26]. We note 
that the same effective Hamiltonian also describes monolayer 
graphene with staggered sublattice potential [15, 16]. This 
is not surprising, as both systems have the same symmetry 
properties. What distinguishes M0S2 from graphene is the 
strong SOC originated from the metal d-orbitals. The con- 
duction band-edge state is made of d^^ orbitals and remains 
spin-degenerate at K points, whereas the valence band-edge 
state splits. Approximating the SOC by the intra-atomic con- 
tribution L • S,WQ find the total Hamiltonian given by 



H = atirkxCTx + kyCfy) + — - Ar 



(3) 



where 2 A is the spin- splitting at the valence band top caused 
by the SOC and 5^ is the Pauli matrix for spin. The spin-up 
(t) and spin-down (|) components are completely decoupled 
and Sz remains a good quantum number. We emphasize that 
the spin splitting does not depend on the model details; it is 
a general consequence of inversion symmetry breaking, sim- 
ilar to the Dresselhaus spin splitting in zinc-blende semicon- 
ductors [27]. Time-reversal symmetry requires that the spin 
splitting at different valleys must be opposite (Fig. Ic) [28]. 

The valley Hall and spin Hall effects are driven by the 
Berry phase associated with the Bloch electrons. It has been 
well established that in the presence of an in-plane electric 
field, an electron will acquire an anomalous velocity propor- 
tional to the Berry curvature in the transverse direction [29], 
giving rise to an intrinsic contribution to the Hall conductiv- 
ity [30], cr^^^ = {e^/h) J[dk]f{k)n{k), where f{k) is the 
Fermi-Dirac distribution function, and [dk] is a shorthand 
for dk/ (2ti)'^. The Berry curvature is defined by Vtnik) = 
z ■ Vfc X {un{k)\i\/ k\un{k)) , where \un{k)) is the periodic 
part of the Bloch function and n is the band index. For mas- 
sive Dirac fermions described by the effective Hamiltonian in 
Eq. (3), the Berry curvature in the conduction band is [15]: 



n^k) 



(4) 



3 



Note that the Berry curvatures have opposite sign in opposite 
valleys. In the valence band, we have: ^v{k) = —^c{k). 
In the same valley, the Berry curvature is dependent on spin 
through the spin-dependent band gap: A' = A — rs^X. The 
curvature is nearly constant in the neighborhood of K points 
since /S.^ atk (Table I). The valley Hall conductivity (in unit 
oie/h) is then: 

and the spin Hall conductivity (in unit of e/2) is: 

a,^ = [dk]{fn^^{k)^ln^^{k) - fn,l{k)Qn,i{k)) , 

where the integration is performed over the neighborhood of 
one K point. For moderate hole doping with Fermi energy 
lying between the two split valence band tops (illustrated by 
the dot-dashed line in Fig. 2b), the valley and spin Hall con- 
ductivities of holes are the same, given by 
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for fi <^ A — A, where /i is the Fermi energy measured from 
the valence band maximum. If the system is electron doped, 
we must consider both conduction bands which are degenerate 
at K points and have small spin- splitting quadratic in k (see 
Fig. 2b). We find that 
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where /i is the Fermi energy measured from the conduction 
band minimum. The spin Hall conductivity is about A/ A of 
the valley Hall conductivity. 

The robustness of the valley and spin Hall effects is closely 
related to the relaxation time of the valley and spin index. 
Flipping of valley index require atomic scale scatters, since 
the two valleys are separated by a wave vector comparable 
with the size of Brillouin zone. Spin flips requires the cou- 
pling with magnetic defects, since 5^ is a good quantum num- 
ber at the conduction and valence band edges. In the conduc- 
tion band, valley scattering could be slow in the bulk at the 
clean limit, but will be facilitated on the boundaries by valley 
mixing except with perfect zigzag edge. In the valence band, 
by the relatively large valley-contrasting spin splitting (~ 0.1 - 
0.5 eV), valley and spin can only be simultaneously flipped to 
conserve energy which require atomic scale magnetic scatters. 
In the absence of such scatters, we expect holes have long spin 
and valley lifetimes both in the bulk and on the boundary. 

Next we look at optical inter-band transitions from the 
spin-split valence band tops to the conduction band bottoms. 
The coupling strength with optical fields of a± circular po- 
larization is given by V±{k) = Vx{k) ± iVy{k), where 
(k) = mo {uc{k) I ^ \uy {k)) is the interband matrix el- 
ement of the canonical momentum operator and mo is the free 
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FIG. 2. (color online). Coupled spin and valley physics in monolayer 
group- VI dichalcogenides. The electrons and holes in valley K are 
denoted by white ' — ', '+' symbol in dark circles and their counter- 
parts in valley —K are denoted by inverse color, (a) Valley and spin 
Hall effects in electron and hole doped systems (see text), (b) Val- 
ley and spin optical transition selection rules. Solid (dashed) curves 
denote bands with spin down (up) quantized along the out-of-plane 
direction. The splitting in the conduction band is exaggerated, uju 
and (jJd are respectively the transition frequencies from the two split 
valence band tops to the conduction band bottom, (c) Spin and val- 
ley Hall effects of electrons and holes excited by linearly polarized 
optical field with frequency Uu. (d) Spin and valley Hall effects of 
electrons and holes excited by two-color optical fields with frequen- 
cies (jJu and Ud and opposite circular polarizations. 



electron mass [16, 31]. For transitions near K points, we find 

A' 
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Since A' ^ atk, the interband transitions are then coupled 
exclusively with a-\- (a—) circularly polarized optical field at 
the K (—K) valley. Optical field couples only to the orbital 
part of the wave function and spin is conserved in the opti- 
cal transitions. By the valley-contrasting spin splitting of the 
valence band tops, the valley optical selection rule becomes 
spin-dependent selection rules, as illustrated in Fig. 2b. Ud 
and Uu denote here the two band edge excitonic transition fre- 
quencies from the spin-split valence band tops (see Fig. 2b). 
Because of the spin- valley coupling and the valley optical se- 
lection rule from inversion symmetry breaking, spin and light 
polarization are related in the opposite ways at the two fre- 
quencies, similar to the interband transition involving heavy 
hole and light hole in III-V semiconductors. One may expect 
a sign reversal for magneto-optical effects such as Faraday ro- 
tation and Kerr rotation by spin-polarized electron when the 
frequency changes from Uu to ujd- 

Selective excitation of carriers with various combination of 
valley and spin index becomes possible using optical fields of 
different circular polarizations and frequencies. Optical field 
with (7+ circular polarization and frequency uJu (uJd) can gen- 
erate spin up (down) electrons and spin down (up) holes in 
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TABLE IL Photo-induced spin, valley and charge Hall effects, cfs 
and (Jv^^^ denote the spin and valley Hall conductivity respectively 
contributed by the photo-excited electrons (holes), a is the total 
charge Hall conductivity from both carriers. All conductivities are 
normalized by the photo-excited carrier density of electron or hole, 
and only intrinsic contribution is considered. 
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valley K, while the excitation in the —K valley is simply 
the time reversal of the above [32]. Such a spin and valley 
dependent selection rule can be used to generate long lived 
spin and valley accumulations on sample boundaries in a Hall 
bar geometry. Consider the photo-excitation of electrons and 
holes, which are then dissociated by an in-plane electric field, 
driving a longitudinal charge current (Fig. 2c-d). The photo- 
excited electrons and holes will also acquire opposite trans- 
verse velocities because of the Berry curvatures in the conduc- 
tion and valence band, and moved to the two opposite bound- 
aries of the sample. This leads to Hall current of valleys, spins 
or charges, depending on the polarization and frequency of the 
optical field. In Table II, we give the signs and order of mag- 
nitude estimation of the valley, spin and charge Hall currents 
in the clean limit. 

Excitation with circular polarizations will generate a charge 
Hall current which can be detected as a voltage. The sign of 
the voltage is exclusively determined by the circular polariza- 
tion and is independent of the frequency. Excitation with lin- 
ear polarizations has more interesting consequences. For ex- 
ample, by excitation with linearly polarized optical field with 
frequency uju, there is a spin Hall current and a valley Hall 
current in the absence of a charge Hall current. Spin up elec- 
trons from the K valley and spin up holes from the —K valley 
are accumulated on one boundary, while their time reversals 
are accumulated on the other boundary (Fig. 2c). Thus, each 
boundary can remain charge neutral while carrying a net spin 
polarization as well as a net valley polarization. Recombina- 
tion of these excess electrons and holes are forbidden by the 
optical transition selection rules unless assisted by processes 
which flip both the valley and spin index. Holes are expected 
to have much longer spin and valley lifetimes on the bound- 
ary. Thus electrons will get unpolarized first and recombine 
with the spin and valley polarized holes, accompanied by the 
emission of photons with opposite circular polarizations on 
the two boundaries. If there is strong valley mixing for elec- 
trons on the boundary, the decay of the overall spin and valley 
polarization is determined by the spin relaxation time of the 
electrons. 

Another interesting excitation scenario is by a non- 



degenerate optical excitation, consisted of a cr+ polarized 
component with frequency 00 u and a a— polarized component 
with frequency 00 d. This will excite spin up electrons and spin 
down holes in both valleys. The spin Hall and charge Hall cur- 
rents from the electrons will largely cancel with those from the 
holes, while the valley Hall currents from electrons and holes 
add constructively (see Table II). The electrons and holes ac- 
cumulated on the same boundary are of opposite spin and val- 
ley indice. When electrons get valley unpolarized, they can 
recombine with the spin and valley polarized holes, accompa- 
nied by the photon emission with polarization and frequency 
{a-\-^uju) on one boundary and {a—^Ud) on the other. This 
process may provide a direct measurement on the valley life- 
time of electrons on the boundary. 

In summary, we have predicted the valley dependent op- 
tical selection rules for interband transitions near i^-points 
in monolayer M0S2 and other group VI transition metal 
dichalcogenides. The spin-orbit interaction from the metal 
d-orbitals further leads to strong coupling of spin and val- 
ley degrees of freedom, which makes possible selective pho- 
toexcitation of carriers with various combination of valley and 
spin indices. We have also predicted the coexistence of valley 
Hall and spin Hall effects in n-doped and p-doped systems, 
and proposed photo-induced spin Hall and valley Hall effects 
for generating spin and valley accumulations on edges. The 
strong spin- valley coupling can further protect each index: 
with the valley dependent spin splitting of 0(0.1 ) eV at the 
valence band top, flip of spin and valley alone is energetically 
forbidden. These effects suggest the potential of integrated 
spintronic and valley tronic applications. In hybrid systems of 
these monolayers with other spintronics materials, spin index 
may be used as a universal information carrier across differ- 
ent materials, while valley index provides a unique ancillary 
information carrier in the monolayers with logic operations 
between the two enabled by the spin- valley coupling. 
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Note added. - Recently, experimental evidences on the op- 
tical selection rules for inter-band transitions at K points are 
reported in monolayer M0S2 [33]. We also note an indepen- 
dent theoretical work discussing the circular dichroism in the 
entire Brillouin zone [34] . 
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